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Wave  Propagation  at  Regional  Distances,  Contract  #  F19628-89-K-0013 


Report  Summary 

Task  Objectives 

Characterization  of  local  and  regional  wave  propagation  in  the  crust  and 
upper  mantle  around  the  three  short-period  Fennoscandian  arrays, 
NORESS.  ARCESS  and  FINESA. 

Technical  Problem 

Short-period,  small-aperture  arrays  are  key  to  reliable  regional  nuclear  test 
monitoring.  The  array  geometries  enable  phase  detection  at  small 
amplitudes  and  allow  determination  of  backazimuth  and  phase  velocity. 
Assuming  a  velocity  model,  event  distance  can  then  be  deduced  from 
relative  arrival  times  and  phase  velocities  of  major  phases.  In  order  to 
obtain  accurate  source  distances  and  avoid  misidentifying  phases  or  missing 
small  arrivals  in  events,  it  is  therefore  necessary  to  know  the  local  and 
regional  velocity  structure  and  its  effect  on  wave  propagation  near  the 
arrays. 

"  Previpus  work  on  events  at  local  distances  from  the  NORESS  array  has 
shown  diat  mislocations  occur  due  to  missed  small  amplitude  Pn  arrivals  and 
mi^id^tiBed  onset  of  Lg.  Furthermore,  earthquakes  at  local  distances  from 
NORESS  do  net  produce  depth  phases  of  significant  amplitudes,  making 
-source-depth  discrimination  dependent  on  the  existence  of  accurate  travel¬ 
time  CiliVes  for  near-surface  sources.  The  present  work  focuses  on 
.,,,^studyin^  regional  wave  propagation  near  ARCESS,  and  in  addition  the 
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effects  of  receiver  structure  at  NORESS  and  ARCESS.  Behavior  of  the 
main  regional  phases  with  distance  and  azimuth  is  studied  in  ARCESS  data 
from  regional  events.  Receiver  structure  under  NORESS  is  studied  through 
receiver  function  analysis  of  teleseismic  events  recorded  at  the  three- 
component  intermediate-period  elements,  NREO. 

General  Methodology 

Regional  phases  are  studied  through  array-analysis,  f-k  analysis  is  applied 
in  sliding  time  windows  to  detect  phases  and  determine  phase  velocities  and 
azimuth  of  approach.  Composite  seismograms  for  the  events  are  made 
from  time  pieces  of  array  beams,  with  each  time  section  representing.an 
arrival.  Composite-seismogram  record  sections  are  constructed  and 
compared  to  synthetic-seismogram  record  sections,  calculated  by 
wavenumber  integration  in  plane-layered  constant-velocity  models. 

Receiver  structure  is  studied  through  receiver  function  analysis  of 
teleseismic  events.  Receiver  functions  are  calculated  for  all  events  and 
stacked  to  lower  the  effects  of  noise.  Synthetic  receiver  functions  are  then 
calculated  for  different  velocity  models  and  matched,  by  trial  and  error,  to 
the  stack  of  observed  receiver  functions. 

Technical  Results 

Mislocations  with  the  Intelligent  Monitoring  System  of  small  events  around 
ARCESS  are  mainly  due  to  misidentified  onset  of  Lg\  multipathing  also 
causes  some  azimuthal  uncertainties.  Mislocations  due  to  missed  Pn  arrivals 
are  not  observed,  since  Pn  amplitudes  are  relatively  large.  Lateral  velocity 
variations  appear  to  be  associated  with  the  Lapland  Granulite  region  and  the 
Polmak-Pasvik-Pechenga  belt,  causing  apparent  multipathing  and  phase 
velocity  variations.  Crustal  thickness  around  ARCESS  varies  between  42 
and  45  km. 

Characteristic  of  the  seismic  record  section  from  ARCESS  is  the  absence  of 
multiple  crustal  arrivals  at  any  distance.  This  can  be  explained  by  a  velocity 
gradient  above  the  Moho,  causing  the  energy  of  each  reverberation  and 
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multiples  to  be  concentrated  over  a  small  distance  interval.  Rg  is  observed 
out  to  at  least  400  km  distance. 

Crustal  thickness  under  NORESS  is  35  km  and  may  be  shallower  towards 
the  Oslo  Graben.  No  large  crustal  discontinuities  are  detected  with  receiver 
function  analysis.  Apparent  secondary  Rg  waves  seem  to  stem  from  the 
surface  Topography  at  Lake  Mjpsa. 

Important  Findings  and  Conclusions 

Significant  lateral  velocity  variations  affect  the  wave  propagation  in  the 
ARCESS  region.  Higher  velocity  gradients  above  and  below  the  Moho  and 
a  thicker  crust  also  result  in  different  record  sections  at  ARCESS  from 
those  at  NORESS:  At  NORESS,  a  low  velocity  gradient  above  the  Moho 
results  in  multiple  crustal  reverberations  observed  at  any  distance,  and  a 
low  velocity  gradient  below  the  Moho  causes  small  Pn  waves.  At  ARCESS, 
on  the  other  hand,  the  high  gradients  produce  larger  Pn  waves  and  focus  the 
energy  of  each  crustal  reverberation  over  a  small  distance  interval.  The 
low  relief  in  the  ARCESS  region  also  allows  more  efficient  propagation  of 
Rg. 

Due  to  the  discrete  nature  of  arrivals  in  the  array  seismograms,  the 
characteristics  of  the  record  section  at  ARCESS,  as  well  as  at  NORESS  can 
be  explained  in  terms  of  plane-layered  velocity  models.  Although  actual 
velocity  and  thickness  variations  cause  the  observed  record  sections  to 
deviate  from  these  simple  models,  the  models  enable  reliable  identification 
of  the  various  phases  observed  and  thus  result  in  more  accurate  locations. 

Significant  Hardware  Development 
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Introduction 


Short-period  seismograms  from  events  at  local  and  regional 
distances  consist  of  a  multitude  of  complex  arrivals  mixed  in  with 
incoherent  scattered  energy  and  noise.  For  source  distances  well  into  the 
regional  range  however,  the  seismograms  are  commonly  dominated  by 
discrete  arrivals  with  distinct  phase  velocities.  Timing,  amplitude  and 
phase  velocities  of  these  dominant  phases  are  controlled  by  crustal  and 
upper  mantle  structure  as  well  as  source  parameters.  In  order  to  properly 
identify  the  arrivals  and  thus  obtain  reliable  locations  for  the  events, 
average  structure  along  the  propagation  paths  must  first  be  obtained. 
When  amplitude  and  travel-time  behavior  of  the  major  phases  with  distance 
is  known  for  all  azimuth  directions,  it  becomes  possible  to  distinguish  near¬ 
source  surface  reflections,  or  depth  phases,  from  the  train  of  arrivals  in  a 
regional  seismogram.  Depth  phases  reveal  source  depth,  which  for 
discriminatory  purposes  is  the  most  important  source  parameter. 

Array  data  allows  the  association  of  phase  velocities  with  individual 
arrivals.  This  makes  phase  identification  more  reliable  and  enables  tracing 
arrivals  over  distance.  From  the  travel-time  curves  obtained,  the  gross 
structure  may  then  be  deduced.  Due  to  the  discrete  nature  of  the  arrivals 
and  the  domination  of  coherent  phases  over  incoherent  scattered  energy,  a 
central-element  seismogram  can  in  most  cases  be  approximated  by  a 
composite  of  array  beams,  with  each  time  section  of  a  beam  representing 
an  arrival.  This  has  been  demonstrated  for  several  events  recorded  at 
NORESS  (Vogfjord  and  Langston,  1991).  The  composite  seismogram 
enhances  coherent  arrivals,  while  incoherent  energy  and  noise  are 
diminished.  The  difference  between  a  composite  seismogram  and  the 
original  data  gives  an  estimate  of  the  incoherent,  scattered  energy  caused 
by  structural  heterogeneities  and  determines  the  appropriateness  of 
modeling  with  plane-layered,  isotropic  velocity  structure. 

Results  from  the  short-period  regional  data  can  be  augmented  by 
receiver  function  analysis  of  teleseismic  events  recorded  at  the  three- 
component  elements  of  the  arrays.  In  particular,  the  velocity  function 
under  a  recording  site  can  be  studied  independent  of  path  propagation, 
allowing  velocity  contrasts  such  as  the  crust-mantle  interface  to  be  located. 
Knowing  the  near-array  structure,  the  site  response  can  then  be  modeled 
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and  separated  from  the  regional  path  effects.  The  receiver-function 
method  entails  deconvolving  the  vertical  component  of  motion  in  the 
teleseismic  P-wave  window  from  the  radial  component.  This  removes  the 
effects  of  source  time  function  and  source  stmcture,  as  well  as  the  effects 
of  instrument  response  and  propagation  through  the  mantle.  What  is  then 
left  on  the  deconvolved  radial  component  are  the  P-to-S  (Ps)  conversions 
occurring  at  velocity  discontinuities  under  the  receiver.  Descriptions  of 
the  receiver-function  method  are  given  by  Langston  (1979)  and  Owens  et 
al.  (1987). 

The  crustal  and  upper  mantle  structure  around  the  NORESS  and 
ARCESS  arrays  has  been  extensively  studied  with  refraction  surveys.  The 
structure  around  NORESS  is  defined  by  three  geological  provinces:  The 
Precambrian  crust  of  the  southern  Baltic  Shield,  the  Oslo  Graben  within 
the  Precambrian  region  and  the  Caledonian  fold  chain,  which  extends  along 
the  entire  west  coast  of  Scandinavia  (Figure  1).  Velocity  models  for  the 
Graben  and  the  Precambrian  region  east  of  the  Graben  were  obtained  from 
a  study  of  local  events  recorded  at  NORSAR  (Gundem,  1984).  Several 
different  velocity  models  have  been  obtained  for  the  Caledonian  region 
(Kanestr0m  and  Haugland,  1971;  SellevoII  and  Warrick,  1971; 
Mykkeltveit,  1980).  The  Caledonian  model  used  here  is  a  combination  of 
an  upper  crust  from  the  Precambrian  model  and  a  lower  crust  similar  to 
the  models  obtained  from  the  Arsund-Otta  (Mykkeltveit,  1980)  and  the  3-4 
(Kanestr0m  and  Haugland,  1971)  profiles  (Figure  1).  NORESS  is  located 
in  the  Precambrian  region  and  therefore  the  Precambrian  model  is  used  in 
the  receiver  function  modeling.  The  local  and  regional  events  presented 
from  NORESS  are  all,  except  one,  located  in  the  Caledonides  with 
propagation  paths  mostly  within  the  region.  The  travel -time  curves  and 
synthetic  seismograms  calculated  for  the  NORESS  region  are  therefore 
based  on  the  Caledonian  model.  The  two  models  are  shown  in  Figure  2a 
and  2b. 

ARCESS  is  located  in  the  oldest,  Archaean,  part  of  the  Baltic  Shield 
(Figure  3),  just  west  of  the  Proterozoic,  Lapland  Granulite  Belt.  The 
granulite  belt  is  a  wedge-shaped  body,  believed  to  be  a  slice  of  continental 
crust  overthrust  on  to  the  Archaean  basement  to  the  south  in  a  continent- 
continent  collision  during  the  Proterozoic.  The  suture  zone  between  these 
colliding  continents  is  along  the  Polmak-Pasvik-Pechenga  Belt,  which 
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consists  of  Proterozoic  volcanic  rocks  and  runs  approximately  NW-SE 
along  the  Kola  peninsula  (Gaal  et  al,  1989).  The  PPP  Belt  has  been 
penetrated  by  the  Kola  superdeep  drillhole  (Figure  3).  The  velocity 
profile  in  the  drillhole  is  reversed,  with  highest  velocities  in  the  upper  6.8 
km,  which  constitute  the  Proterozoic  section.  Below  that  the  Archaean 
basement  is  reached  (Kozlovsky,  1987). 

Several  refraction  profiles  cross  the  region:  The  northernmost 
section  of  the  Fennolora  profile  runs  just  west  of  ARCESS  (Galson  and 
Mueller,  1986;  Guggisberg  et  al.,  1991;  Lund,  1987)  and  the  Polar  profile, 
crosses  the  Granulite  belt  and  the  suture  zone  east  of  ARCESS  (Luosto  et 
al.,  1989)  (Figure  3).  The  Polar  profile  also  has  a  reflection  section  where 
the  profile  crosses  the  Lapland  Granulite  Belt  (Behrens  et  al.,  1989).  Kola 
Peninsula  has  been  explored  on  several  refraction  profiles  (Azbel  et  al., 
1989;  Glaznev  et  al.,  1989)  and  the  area  south  of  the  array  is  sampled  by 
the  Finlap  profile  (Luosto  et  al,  1983).  One-  (Lund,  1987)  and  two- 
dimensional  (Guggisberg  et  al.,  1991)  velocity  interpretation  of  the 
northern  section  of  Fennolora  includes  low-velocity  zones  in  the  upper 
crust  of  the  Caledonian  region  and  a  45  km  thick  crust .  2-D  modeling  of 
the  Polar  profile  reveals  a  NW  dipping  high-velocity  layer  coinciding  with 
the  Granulite  belt,  and  an  undulating  Moho,  with  the  thinnest  crust  in  the 
middle  of  the  profile.  2-D  interpretation  of  the  Kola  Peninsula  profiles 
requires  an  approximately  40  km  wide  and  10  km  deep  high-velocity 
gradient  zone  in  the  approximate  location  of  the  Granulite  belt  and  an 
undulating  Moho  (Azbel  et  al.,  1989).  Moho  depths  increase  from  36  km, 
under  northeastern  Kola  Peninsula,  to  44  km  under  the  southwestern 
Peninsula.  One-dimensional  modeling  of  Finlap  includes  a  50  km  thick 
crust  (Luosto  et  al.,  1983). 

Location  :  of  the  events  selected  are  shown  in  Figure  3.  All  events 
west  of  ARCESS  are  within  the  Archaean  region,  with  partial  travel-paths 
through  the  Granulite  belt,  and  most  events  on  the  Kola  peninsula  are 
located  near  the  PPP  zone.  Other  events  are  in  the  younger,  Proterozoic 
region  south  of  the  array,  and  one  is  in  the  Caledonian  region,  which 
extends  along  the  western  coast.  The  velocity  model  used  here,  in  travel¬ 
time  and  synthetic  seismogram  calculations  is  a  plane-layered 
approximation  to  the  northern  section  of  the  Polar  profile  and  is  shown  in 
Figure  2c. 
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In  a  previous  study  local  and  regional  events  recorded  at  NORESS 
were  analyzed  (Vogfjdrd,  1991).  In  the  present  study  the  research  is 
focussed  on  regional  ARCESS  data.  Arrivals  in  the  seismograms  are 
detected  and  identified  with  f-k  analysis  in  sliding  time  windows,  beams  are 
formed  for  the  major  phases,  and  time  sections  of  the  beams  then  pieced 
together  to  form  the  composite  seismogram  for  each  event.  A  seismic 
record  sections  of  composite  seismograms  from  ARCESS  is  constructed 
and  compared  to  a  record  section  from  NORESS,  and  behavior  of  the 
dominating  phases  with  distance  is  examined.  The  composite  record 
sections  reveal  the  behavior  of  relative  amplitudes  and  travel  times  with 
distance,  thus  allowing  conclusions  about  crustal  and  upper  mantle 
structure  to  be  drawn.  Synthetic  P-wave  record  sections  calculated  for 
plane-layered  isotropic  models  ire  also  presented  in  an  attempt  to  explain 
the  apparent  differences  between  the  two  regions. 

Receiver  functions  from  teleseismic  events  recorded  on  the  three - 
component,  intermediate-period  chaimels  EO,  at  the  center  of  NORESS  and 
ARCESS  are  being  analyzed.  Results  from  NORESS  are  presented  and  the 
Ps  conversion  at  the  Moho  is  modeled.  Crustal  discontinuities  are  not 
resolved,  but  a  surface  scatterer  is  identified,  producing  strong  coda  waves. 


Research  Accomplished: 

Regional  Pnases  at  ARCESS 

Twenty  two  seismic  events  at  local  and  regional  distances  from 
ARCESS  were  selected  from  the  IMS  bulletin  at  CSS.  The  events  cover  a 
distance  range  of  175  to  575  km  and  local  magnitude  range  (Af;)  is  from 
2.0  to  3.2.  Most  locations  coincide  with  known  mines  and  are  therefore 
likely  mine  explosions.  A  study  of  eight  events  has  already  been 
compleded.  Event  locations  are  shown  in  Figure  3.  Due  to  misidentified 
onset  of  Lg.two  events  are  mislocated  by  the  Intelligent  Monitoring 
System,  by  approximately  30  km.  In  event  167198  at  175  km  distance,  Sg 
is  the  first  arrival  followed  by  a  larger  SmS  phase  (Figure  4),  which  is 
mistakenly  picked  as  the  first  S-arrival.  In  event  335253  at  220  km 
distance  SmS  is  missed  and  the  following  larger  phase,  2xSmS,  is  picked  as 
the  onset  of  Lg.  Events  due  east  from  ARCESS  show  a  significant  amount 
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of  variation  in  backazimuth  within  each  phase.  This  variation,  which  is 
probably  due  to  lateral  multipathing  may  be  causing  a  mislocation  of  up  to 
10’  in  event  326360  at  200  km  distance. 

Composite-seismogram  record  sections  of  the  events  completed  are 
shown  in  Figure  4a  and  b,  with  stacking  velocities  indicated  above  each 
beam  section.  Event  348962  has  been  high-pass  filtered  at  1.5  Hz  to  reduce 
low-frequency  noise.  The  travel-time  curves  shown,  were  calculated  for  a 
plane-layered  approximation  to  the  Polar-Profile  model  (Luosto  et  al., 
1989)  with  a  crustal  thickness  of  42  km.  The  move-out  of  the  dominant 
phases  approximately  follow  the  travel-time  curves,  indicating  an  adequate 
velocity  model.  Lack  of  events  within  175  km  from  the  array,  however, 
precludes  resolution  of  detailed  crustal  structure  and  the  observed 
deviations  of  the  Moho  reflections  from  the  travel-time  curves  also 
indicates  variations  in  crustal  thickness. 

The  main  characteristics  of  the  record  section  in  Figure  4a  are 
apparent  discrete  arrivals,  each  mostly  dominated  by  a  single  phase 
velocity.  In  event  240323,  at  575  km  distance,  the  phase  velocities 
obtained  in  the  Lg  wave  train  however  are  more  varied,  although  the  wave 
train  is  dominated  by  3.0  km/s  and  4.6  km/s  waves.  Relatively  large  upper 
mantle  waves,  Pn  and  Sn,  are  also  observed,  probably  caused  by  a 
considerable  velocity  gradient  just  below  the  Moho.  Pg  is  the  first  arrival 
in  event  167198  at  175  km  distance,  at  200  km  Pg,  PmP  and  Pn  arrive 
approximately  simultaneously,  and  beyond  that  Pn  is  the  first  arrival.  The 
relative  arrival  times  of  the  inferred  Pg  and  PmP  phases  in  the  event  at  175 
km  distance,  as  compared  to  the  travel-time  curves,  may  indicate  a  deeper 
Moho  at  90  km  distance  east  of  the  array.  However,  it  could  also  be  caused 
by  a  ripple  fired  explosion,  in  which  case  both  P  arrivals  are  interfering  Pg 
and  PmP.  Event  348962  also  appears  to  be  near  the  crossover  distance  for 
Pn,  PmP  and  Pg  indicating  a  deeper  Moho  than  42  km. 

The  first  Moho  reflections  (PmPISmS)  dominate  the  seismograms  out 
to  approximately  300  km  distance.  Beyond  that  distance  the  second  crustal 
reverberations  take  over,  while  the  first  ones  disappear.  The  upper  mantle 
waves  experience  increase  of  amplitude  with  distance  so  that  in  event 
240323,  at  575  km  distance  they  dominate.  A  clear  double  arrival  is 
observed  in  the  Pn  wave  train  at  575  km  distance  followed  by  a  third,  9.1 
km/s  arrival,  possibly  2xPn.  Whether  the  second  arrival  is  a  source  or  a 
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structure  effect  will  be  determined  by  examining  other  events  from  this 
location.  Only  in  event  3352S3,  at  350  km  distance  are  there  more  than 
one  Moho  reverberation  observed,  i.e.  SmS  and  2xSmS.  This  absence  of 
multiples  can  be  explained  by  a  strong  velocity  gradient  in  the  lower  crust, 
causing  focussing  of  the  energy  of  each  cmstal  reveiberation  over  a  short 
distance  interval.  This  is  consistent  with  the  Polar-Profile  model  (Luosto 
et  al.,  1989). 

Rg  is  observed  out  to  400  km  distance,  and  moves  out  with  a  group 
velocity  of  2.9  km/s.  The  interference  apparent  in  Rg  of  events  167198  at 
86.0°  backazimuth  and  event  326360  at  95.5°  backazimuth,  is  not  observed 
in  Rg  of  event  282554  (Figure  4b),  which  is  at  approximately  the  same 
location  as  event  326320  and  is  also  much  larger  than  the  other  two.  This 
difference  therefore  is  probably  a  source  effect  -  such  as  a  ripple-fired 
explosion  -  rather  than  a  structure  effect,  although  the  location  of  the 
events  in  and  around  the  PPP  Belt  may  cause  multipathing  near  the  surface. 

Event  392876  is  located  in  the  Precambrian  region  southwest  of 
ARCESS.  The  propagation  path  is  along  the  eastern  margin  of  the 
Caledonian  region  and  crosses  a  section  of  f  '^ennolora  profile  which  has 
complex  upper  crustal  structure  (Lund,  1987;  Guggisberg  et  al.,  1991). 
This  may  be  the  cause  of  the  small  arnplitude  Pn  and  the  apparent  multiple 
crustal  arrivals  observed  in  the  P  wave  train.  The  multiple  arrivals 
however,  could  be  due  to  source  time  function  and/or  source  depth.  Rg 
was  not  detected  in  the  unfiltered  records,  but  due  to  the  low  amplitudes, 
may  have  been  buried  in  the  low-frequency  noise.  The  apparent  absence  of 
Rg  could  be  the  result  of  considerable  source  depth,  but  the  location  of  the 
event  near  Kiruna,  suggests  an  explosion  source.  Analysis  of  more  events 
from  this  location  is  necessary  to  distinguish  stmcture  and  source  effects. 

The  record  section  in  Figure  4a  shows  a  distinctly  different 
character  from  record  sections  of  composite  seismograms  from  the 
NORESS  array.  For  comparison,  a  record  section  of  seven  events  from  a 
previous  study  (Vogfjbrd,  1991),  with  locations  as  shown  in  Figure  1  is 
plotted  in  Figure  5.  Travel-time  curves  calculated  for  a  Caledonian  model, 
36-km  thick,  are  superimposed  on  the  record  section.  All  events,  except 
one,  B-5,  are  located  within  the  Caledonides,  with  propagation  paths 
mainly  within  Caledonian  structure.  The  record  section  reveals  small 
upper  mantle  waves,  more  than  one  Moho  reverberation  in  events  from 
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distances  beyond  approximately  200  km,  and  inefficient  propagation  of  Rg 
waves  within  the  Caledonides.  The  small  upper-mantle  waves  are  caused 
by  a  small  or  no  velocity  gradient  below  the  Moho.  The  occurrence  of 
large  double  reverberations  (2xPmP  and  2xSmS),  while  the  first  reflections 
(PmP  and  SmS)  still  remain  large  can  be  explained  by  gradual  velocity 
increase  in  the  upper  to  intermediate  crust  and  a  near-constant  velocity  in 
the  lower  crust.  Due  to  scattering  by  surface  topography  in  the 
Caledonides,  Rg  does  not  propagate  to  200  km  distance.  Propagation 
through  Precambrian  structure  is  more  efficient,  as  can  be  seen  in  event  B- 
5,  which  has  an  Rg  wave  of  significant  amplitude.  The  model  used  for 
travel-time  calculations  represents  average  structure  in  the  Caledonides:  a 
plane-layered  36  km  thick  crust,  with  a  gradual  velocity  increase  in  the 
intermediate  crust  and  constant  velocity  in  the  lower  crust  (Vogfjdrd, 
1991). 

To  illustrate  the  difference  between  the  Caledonian  and  the  Polar- 
profile  velocity  models,  synthetic  seismograms  are  calculated  by 
wavenumber  integration  in  plane-layered  constant-velocity  models,  which 
are  approximations  to  the  gradient  models.  The  P-wave  velocity  models 
are  plotted  in  Figure  2b  and  c.  To  facilitate  the  comparison,  the  thickness 
of  the  bottom  crustal  layer  in  the  Caledonian  model  is  increased  to  make 
crustal  thickness  equal  that  of  the  Polar-profile  model.  Integration  is 
performed  over  ray  parameters  from  0.01  to  0.2  s/km;  including  only  P- 
waves.  The  source  time  function  used  is  triangular,  of  0.1  s  duration  and 
sampling  frequency  is  10  Hz.  The  synthetic  record  sections  together  with 
travel-time  curves  are  plotted  in  Figures  6  and  7.  The  synthetic  Polar- 
Profile  record  section,  in  Figure  6,  shows  large  amplitude  Pn  waves  and 
concentration  of  the  energy  in  each  Moho  reverberation  over  a  distance 
interval  of  roughly  200  km.  This  is  in  general  accordance  with  the 
composite  record  section  in  Figure  3a.  The  Caledonian-model  record 
section  in  Figure  7,  on  the  other  hand,  has  small  Pn  waves  and  two  to  three 
Moho  reverberations  of  significant  amplitude  at  each  distance  beyond  the 
critical  distance  for  2xPmP.  Due  to  the  gradual  increase  in  velocity 
throughout  the  intermediate  crust,  only  the  Moho  reflections  show  up  as 
discrete  impulsive  arrivals. 
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Receiver  Functions  at  NORESS 
21  teleseismic  events  recorded  at  NORESS  were  used  to  calculate 
receiver  functions  in  order  to  obtain  crustal  thickness  and  velocity  profile 
under  the  array.  The  events,  listed  in  Table  1  are  at  epicentral  distances  of 
30“  to  90",  with  Mb  >  5.8.  Azimuthal  coverage  is  mostly  in  the  northern 
quadrants.  Event  locations  are  shown  in  Figure  8.  Due  to  noise  in  the  data 
stacking  was  necessary  to  improve  the  signal-to-noise  ratio.  In  the  stacked 
data,  the  only  observed  Ps  conversion  occurs  at  the  Moho.  Using  the 
Precambrian  velocity  model,  synthetic  receiver  functions  were  calrulated 
and  fitted  to  the  stacked  receiver  function.  By  trial  and  error,  the  best  fit 
was  obtained  for  a  35  km  thick  crust.  In  accordance  with  the  stacked  data, 
the  synthetic  receiver  function  produced  by  the  Precambrian  model  has 
only  one  apparent  Ps  conversion,  which  occurs  at  the  Moho.  The 
discontinuities  in  the  crust  are  too  small  to  produce  significant  Ps 
conversions.  The  data  and  synthetic  are  shown  in  Figure  9a.  Two 
individual  events  were  of  high  enough  quality  to  allow  inspection  of  Moho 
depth  variation  with  backazimuth.  TTie  receiver  functions  for  these  events 
and  synthetic  are  shown  in  Figure  9b  and  c.  For  a  backazimuth  of  336*, 
the  Moho  depth  is  inferred  35  km,  but  for  a  backazimuth  of  262’,  a  Moho 
depth  of  33  km  is  deduced.  Variations  in  crustal  thickness  around  the 
NORESS  array  were  indicated  by  previous  studies  (Berteussen,  1977, 
Gundem,  1984),  with  thickness  decreasing  under  the  Oslo  Graben.  This  is 
consistent  with  the  receiver-function  results.  Event  22,  from  75’  distance 
and  38’  backazimuth  has  a  high  s/n  ratio  and  large  coda  waves  trailing  the 
first  P  arrival  by  10  seconds  (Figure  10a).  The  particle  motion  in  the 
propagation  plane  from  the  event  is  shown  in  Figure  10b  for  the  15-20 
second  time  interval.  If  the  horizontal  components  are  rotated  to  45’  and 
135’,  as  shown  in  Figure  11a,  the  coda  wave  is  confined  to  the  vertical  and 
45’  components.  The  particle  motion  in  the  vertical-45’  plane,  shown  in 
Figure  11b,  is  retrograde  elliptical.  The  wave  is  therefore  a  secondary  Rg 
wave  caused  by  the  surface  topography,  approximately  27  km  SW  of 
NORESS.  This  location  coincides  with  lake  Mjpsa. 
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Conclusions  and  Future  Plans: 


Location  of  the  ARCESS  array  just  east  of  the  Lapland  Granulite 
belt  appears  to  have  an  effect  on  the  regional  waveforms  observed. 
Variations  in  backazimuth  are  considerable  for  events  located  near  the 
Pechenga  region  and  phase  velocities  in  the  Lg  wavetrain  from  156° 
backazimuth  show  considerable  variation. 

A  significant  difference  is  observed  between  the  short-period 
seismograms  recorded  at  NORESS  and  ARCESS,  from  local  and  regional 
events.  The  difference  can  be  explained  by  a  greater  velocity  gradient  in 
the  lower  crust  and  upper  mantle  of  northern  Fennoscandia  than  in 
southern  Norway.  Furthermore,  due  to  the  lower  relief  in  the  area  south 
and  west  of  ARCESS,  the  Rg  wave  propagates  to  greater  distances  than  in 
southern  Norway.  Mislocations  by  NORESS  occur  due  to  missed  Pn 
arrivals  and  misidentified  onset  of  Lg  (Vogfjdrd  and  Langston,  1991).  At 
ARCESS,  Pn  is  larger  and  therefore  more  easily  detected.  Consequently, 
mislocations  by  ARCESS  are  mainly  due  to  misidentified  Lg  onset.  The 
apparent  multipathing  from  easterly  azimuths  may  also  be  causing  some 
azimuthal  mislocation. 

In  the  continuing  study  coverage  at  local  distances  around  ARCESS 
is  being  improved.  Up  until  recently  only  data  at  distances  >  170  km  and 
of  magnitude  >  2.0  were  saved  for  the  IMS  database  at  CSS.  Data  from 
smaller  distances  are  now  available  on-line  in  the  IMS2  database.  A  scan 
of  the  parameter  data  indicates  that  coverage  at  local  distances  can  be 
significantly  improved.  A  few  good-  to  high-quality  local  events  have 
already  been  retreived  from  IMS2.  A  more  complete  event  coverage  is 
expected  to  allow  resolution  of  the  crustal  and  upper  mantle  velocity 
profile  around  ARCESS.  Modeling  of  Rg  wave  dispersion  should  also  be 
possible  with  an  increased  number  of  local  events. 

Crustal  thickness  under  NORESS  as  obtained  by  receiver-function 
analysis  is  35  km,  with  indications  of  shallower  (33  km)  depths  towards  the 
Oslo  Graben.  The  surface  topography  at  Lake  Mjpsa  scatters  the 
teleseismic  P  waves  producing  a  large  secondary  Rg  wave  arriving 
approximately  10  seconds  later  than  the  P  wave.  The  lake  Mjpsa  scatterer 
has  been  observed  in  previous  teleseismic  studies  (Gupta  et  al.,  1991, 
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Wagner  and  Langston,  1991)  and  apparent  secondary  Rg  waves  have  also 
been  observed  from  a  local  event  (event  D2  in  Figure  1)  with  a 
propagation  path  across  Lake  Mj0sa  (Vogfjbrd,  1991). 

Additional  teleseismic  data  recorded  at  NORESS  and  ARCESS  have 
been  received  from  Norway  and  analysis  of  receiver  functions  at  NORESS 
and  ARCESS  is  now  in  progress. 

The  new  IMS2  database  also  contains  waveform  data  from  the 
FINESA  array  in  Finland.  Analysis  of  local  and  regional  events  recorded 
by  FINESA  is  the  subject  of  future  investigation  in  order  to  compare  wave 
propagation  through  the  Precambrian  crust  in  Southern  Finland  to  the 
Archaean  crust  of  Lapland  and  the  Precambrian  and  Caledonian  crust  of 
Southern  Norway. 
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Relocated  Events 
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Figure  1.  Locations  of  events  recorded  at  NORESS  (triangles).  Events  presented  here 
are  ^own  by  solid  triangles.  The  Caledonian  region  (light)  and  the  Oslo  Graben  (dark) 
are  shaded,  and  locations  of  refraction  profiles  are  shown  by  thick  lines.  Lake  Mjosa  is 
located  approximately  30  km  SE  of  NORESS. 
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Figure  2.  P-wave  velocity  nxxlels:  Models  used  in  travel-time  and  receiver-function 
analysis  (thin  lines)  and  their  constant-velocity  approximations  (thick  lines)  used  in 
synthetic-seismogram  calculations,  a)  Caledonian  model  representing  average 
structure  in  the  Caledonides  near  NORESS;  b)  Precambrian  model  representing 
Precambrian  structure  around  NORESS;  c)  Polar-prorile  model  representing 
average  structure  in  the  Archaean  region  around  ARCESS. 
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Figure  3.  Locations  of  selected  events  (triangles)  recorded  at  ARCESS  (square).  Dashed 
triangles  represent  events  relocated  to  the  positions  indicated  by  the  arrows.  Events 
presented  here  are  shown  by  solid  circles.  Refraction  profiles  are  indicated  and  geologic 
regions  are  shaded:  Caledonian  region  (light-grey);  Lapland  Granulite  Belt  (grey);  Polmak- 
Pasvik-Pechenga  Belt  (dark-grey).  The  cross  in  the  PPP  Belt  marks  the  location  of  the 
super-deep  borehole. 
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Polar  Profile  model  —  42  km  depth  to  Moho 
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Figure  4a  and  b.  Connposite  seismograms  for  the  studied  events  recorded  at  ARCESS,  normalized  to  maximum  amplitude. 
Stacking  velocities  used  to  create  beams  are  indicated  above  each  time  section  in  the  seismo^ms.  Travel-time  curves  of 
major  phases,  calculated  for  the  Polar-profile  model,  are  superimposed  on  the  plots.  Backazimuth  to  the  event  and  (Cont.) 


Polar  Profile  model  -  42  km  depth  to  Moho 


Figure  4  (cont.)  maximum  amplitude  are  shown  on  the  right  side  of  each  trace.  Note  the  apparent  interference  in  the  Rg 
waves  of  events  167 198  and  326360  (record  section  a).  This  is  not  seen  in  Rg  of  event  282554  (record  section  b),  which 
is  much  larger  and  at  approximately  the  same  location  as  event  326360. 
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Caledonian  model  -  42  km  depth  to  Moho 


time  curves  are  superimposed  on  the  plot.  The  synthetics  show  similar  features  as  those  observed  in  the  composite 
seismograms. 


Pole  Profile  model  -  42  km  depth  to  Moho 


letic  P-wave  record  section  for  the  Caledonian  model,  calculated  by  wavenumber  integration.  Travel 
superimposed  on  the  plot.  The  record  section  is  significantly  different  from  Figure  6  and  has  similar 
IS  the  NORESS  composite  record  .section. 


Events  for  Mb>  =  5.a  and  dis  30~90  deg  from  NORESS 
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Receiver  Functions  of  Stacked  Signals  and  H  =  35  km  Precombrion  Model 


Receiver  Functions  of  Jun  11,  1986  event  and  H  =  33  km  Precombrion  Model 


Figiav  9.  Receiver  functions,  for  datt  (thin  line)  and  synthetics  (thick  line),  a)  Receiver 
function  of  stacked  data  and  synthetic  receiver  function  fw  a  35  km  thick  Precambrian 
crust;  b)  Receiver  function  for  event  18  and  the  same  synthetic  as  in  a;  c)  Receiver 
function  for  event  IS  and  synthetic  receiver  function  for  a  33  Ion  thick  PrKambrian  crust. 


Seismograms  (z,  r.  and  t)  for  Apr.  3,  1988  event 


Figure  10  a).  Three-component  plot  (vertical,  radial,  tangential)  of  event  22,  showing  P 
arrival  at  5  sec  and  a  following  coda-wave  on  all  components  at  15  sec. 


Hgure  10  b).  Particle  motion  in  the  propagation  plane  between  15  and  20  sec. 
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Seismograms  (z,  45,  135  deg.)  for  Apr.  3,  1988  event 


Figure  1 1  a).  Three-component  plot  of  event  22,  with  horizontal  components  rotated  to 
45  and  135  deg. 

Particle  Motions  between  15-20  sec 


Figure  1 1  b).  Particle  motion  in  the  veitical-4S  ^g.  plane,  between  15  and  20  sec,  showing 
rettqpade  elliptical  particle  motion. 
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TABLE  1.  Event  Pamneten  for  NORESS  Data 


ID 

Date 

Tune 

UT 

Lat 

deg 

Lon. 

deg 

Af» 

Depth 

km 

Dis. 

deg 

Back  Azi. 

deg 

1 

Mar.  21.  87 

10:50:  35J 

51008 

-177.444 

6.0 

97 

67.4 

6.0 

2 

Feb.  18. 87 

00:00:  52.6 

5151 

-179.31 

61 

33 

68.0 

7.3 

3 

Feb.  7.88 

18:15:  05.6 

50.785 

173.465 

61 

33 

67.9 

12.3 

4 

Feb.  29.  88 

05:31:  41.4 

55.149 

167.430 

6.1 

33 

62.9 

15.3 

5 

Dec.  28.84 

10:46:  135 

56.194 

163.460 

61 

33 

61.3 

17.5 

6 

May.  16.  86 

17:00:  16.7 

47.204 

154.109 

5.8 

25 

681 

26.5 

7 

Jan.  14.  87 

11:03:485 

42.56 

142.88 

65 

99 

69.7 

36.3 

8 

May.  7.87 

03:15:  155 

46.740 

139.230 

6.0 

430 

64.8 

37.0 

9 

Jun.  24.86 

02:53:  09.4 

34.733 

140.502 

6.1 

48 

76.1 

41.3 

10 

Nov.  14.  86 

21:20:  4.7 

23.963 

121.817 

6.1 

33 

78.7 

61.1 

11 

Aug.  26.  86 

09:43:  42 

37.686 

101.412 

6.1 

33 

57.9 

69.4 

12 

Jun.  20.  86 

17:12:  46.5 

31.221 

86.862 

6.0 

33 

561 

85.6 

13 

JuL  7.86 

16:26:  56.9 

10.421 

56.764 

6.4 

10 

601 

126.4 

14 

Mar.  10.  88 

06:17:  235 

10.402 

-60587 

61 

56 

713 

259.4 

IS 

Jua  11. 86 

13:48:  03.3 

10.602 

-61949 

6.0 

34 

731 

261.6 

16 

Jul.  21. 86 

14:42:  265 

37.537 

•118.447 

6.0 

9 

73.9 

320.7 

17 

Mar.  2S.  88 

19:36:  46.4 

61154 

-124.182 

6.1 

10 

519 

335.7 

18 

Oct  5,85 

15:31:  215 

61237 

•124.266 

65 

10 

52.8 

335.8 

19 

Dec.  23. 85 

05:23:  225 

61222 

•124.239 

6.4 

6 

525 

335.8 

20 

Mar.  6.88 

22:35:  38.1 

56.953 

-143.032 

6.8 

10 

61.0 

344.4 

21 

Jun.  19.  86 

09:09:  125 

56.418 

-151729 

6.0 

33 

625 

350.2 

22 

Apr.  3.  88 

01:33;  055 

49.917 

78.945 

6.1 

0 

38.1 

75.1 
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